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DETECTION OF NONLINEAR TRANSFER FUNCTIONS 

BY THE USE OF GAUSSIAN STATISTICS 

By Jack G. Sheppard 
Manned Spacecraft Center 


SUMMARY 


. Failure detection 'will be an essential part of electronic systems in future space 
programs. It would be highly desirable to be able to detect incipient failures while 
equipment is on-line. The possibility of using the statistics of on-line signals as an 
indicator of incipient failures is discussed in this report. As a part of this discussion, 
the concepts of random variables, functions of random variables, and stochastic proc- 
esses are defined in a limited sense. A nonlinearity test that uses ratios of the mo- 
ments of a Gaussian random variable is developed and presented. The results of this 
investigation are encouraging, and the results indicate that further work should be 
pursued. The next logical step would be to apply nonsupervised learning theory to 
determine the statistics of nonstationary input signals and to use the results to detect 
nonlinearities in the output signals of electronic systems. 



INTRODUCTION 



In die manned spacecraft program, great concern exists about operational reli- 
ability; no effort is spared to ensure that, even if the mission cannot be completed 
satisfactorily, the flight crew is returned to earth safely. One of the foundations of 
the operational- reliability philosophy is redundancy. Many kinds of redundancy are 
used (e. g. , redundant testing, redundant inspection, redundant functions, and redun- 
dant equipment). All these techniques are used to ensure that only reliable equipment 
is installed in a spacecraft and that, if any item fails, another piece of equipment or 
mode of operation is available to replace it. 


In the past, the greatest emphasis has been placed on exhaustive testing before 
launch. Usually, the equipment experienced more hours of testing than were experi- 
enced in flight. Degraded operation was permitted in the backup modes, but a signifi- 
cant failure caused immediate mission termination. Because of the nature of future 
manned spacecraft programs, emphasis will be shifted more toward mission- success 
techniques. For space- station operations, permanence will be emphasized, and space 
shuttle vehicles with 100- mission lifetimes will spend more time in space than on the 
ground. Emphasis will be on techniques that minimize ground testing and that avoid 
mission termination because of equipment malfunction. Degraded backup modes of 
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function of t 

first derivative of a with respect to t 
function of t 

first derivative of 0 with respect to t 

i ttl central moment of X 

standard deviation 

second central moment or variance 


RANDOM VARIABLES 

An experiment may be considered in which a marble is selected from a jar con- 
taining marbles of different colors. If a number is assigned to each color, then the 
experiment is as follows. 

1. A marble is selected. 

2. The color of the marble is determined. 

3. The number that corresponds to that color is logged. , 

Thus, an experiment is performed, an outcome is observed, and a number is assigned 
to die outcome. This statement is the definition of a random variable (ref. 1). A ran- 
dom variable is a function from a set of outcomes of an experiment to the set of real 
numbers. 

A different experiment might be the selection of a number from the set of all real 
numbers. In this case, the functional value of the random variable could be the outcome 
of the experiment itself. This is the type of random variable considered in this report. , 
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FUNCTIONS OF RANDOM VARIABLES 


Once a random variable has been defined, the range of the random variable (the 
set of real numbers associated with the experiment outcomes) can be manipulated by 
any one of a multitude of functions. The range may be incremented by a constant, mul- 
tiplied by a constant, subjected to a polynomial transformation, or whatever the imagi- 
nation can devise. Linear transformations, polynomial transformations, and the 
densities and moments of random variables will be discussed in this report. 
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Let x be the value of the random variable X. Then, a linear transformation of 


X is 


. ^{xJsajX+ag 


( 1 ) 
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A polynomial transformation of X is 


g x fc) * a^ + agX + a^ + 


ux 
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n-1 


( 2 ) 


The density of a random variable is a function that describes the relative fre- 
quency of occurrence of the x-values in the range of the random variable. Many types 
of densities exist. The uniform density, for example, states that the relative frequency 
of occurrence of the x-values is constant over a range. This report is concerned with 
the Gaussian or normal density, which is described by the function 
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where a is a spreading factor; Hie graph 
of the Gaussian density function is shown 
in figure 1. - . • •' 

As indicated in figure 1;, the x-values'lpl|ii ,r '•>,*» - 
tend to cluster around a central point, .with 
a decreasing frequency of occurrence as a 
function of distance from this point. 


The probability that X takes on a 
value less than or equal to x^ is called 

the distribution function and is expressed 
by 





Figure J. - Gaussian density function. 



The actual probability that the outcome of the experiment will be between x. and 
x 2 (or equal to either) is found by 1 


P(Xj a x s x 2 ) = 


This integral cannot be evaluated by normal means; therefore, numerical techniques 
must be employed. However, the normal density has been evaluated thoroughly, and 
tables of values are available in almost any text dealing with statistics (ref. 2). 

Moments indicate where the functional values of a random variable are located 

and how the values are spread in the set of real numbers. The equation for the 1 th 
moment of the random variable X is 




( 6 ) 


Central moments are formed by subtracting m, from each value of x in equation (6); 
that is, * 



m^ f^(x) dx 


( 7 ) 


The first moment, called the mean, indicates where the functional values of X are 
centered in the set of real numbers. The second central moment, called the variance, 


o 

is usually designated by the symbol o and Indicates how widely dispersed the func- 
tional values of X are. Higher order moments also indicate dispersion of the 
x values. 


LINEAR TRANSFORMATION OF GAUSSIAN RANDOM VARIABLES 
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where X is a Gaussian random variable. If a is positive, then 


p(y ^yj) = p(j 
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* if equation (10) Is differentiated with respect to y, 
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If a is negative, 
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If equation (IS) is differentiated with respect to y, 
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Thus, by combining equations (11) and (14), 
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Two important points become evident as a result of the preceding discussion. First, 
the linear transformation of a Gaussian random variable results in a Gaussian random 
variable and, second, the variance of the resulting random variable can be expressed 
as the variance of the previous random variable multiplied by the square of the trans- 
formation constant. 


MOMENTS OF A GAUSSIAN RANDOM VARIABLE 
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Leibniz'* rule for differentiation states Oat if 
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By applying Leibniz's rule to equation (18), equation (22) is obtained. 
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If equation (18) is differentiated n times with respect to x, 
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If both sides of wpatioa (26) are divided by then 
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meSods, it should be possible to detect a nonlinear transfer at a Gaussian rai« 
variable. 
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STOCHASTIC PROCESSES 
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A precise definition of stochastic processes has many ramifications. For the 
purposes of this report, a stochastic process will be considered to be a continuous rep- 
etition of the experiment discussed in the section entitled "Random Variables"; that is, 
at all instants of time, the experiment is being performed, and an outcome is available. 
Thus, at any instant of time, a random variable exists, and a range value can be ob- 
tained as a function of the distribution of that random variable. An example of a sto- 
chastic process could be the voltage across the terminals of a battery. If the battery Is 
charged, the voltage has little variation; therefore, the stochastic process is not very 
interesting. A more interesting stochastic process is the voltage across the terminals 
of an antenna. This voltage is the sum of various kinds of noise and the different kinds 
of communications signals currently in use. Hie contrast between these two kinds of 
stochastic processes illustrates an important property of some stochastic processes — 
stationarity. The battery voltage does not change; it remains stationary. The antenna 
voltage moves about drastically. Actually, the process need not sit still to be station- 
ary. A process is strict-sense stationary if its statistics are not affected by a shift in 
the time origin. It is wide- sense stationary if its mean is constant and its autocorrela- 
tion is a function of the time separation coly. Because a Gaussian process is uniquely 
determined by its first two moments, a wide- sense stationary Gaussian process is also 
strict-sense stationary and has a constant mean and variance. This report is primarily 
concerned with Gaussian processes. 


GAUSSIAN PROCESS SAMPLING 


If the normal signal of some linear system is assumed to be a stationary Ga uss ia n 
stochastic process, then, at any point in time, the signal value is the functional value 
of a Gaussian random variable with a fixed mean and variance. Such a random variable 
from a stochastic process is called a sample. By taking many samples from a process, 
inferences can be drawn concerning the distribution of toe source process. The more 
samples drawn, the better toe inferences will be. In this report, toe possibility is 
investigated that very small amounts of distortion in a supposedly linear transfer of a 
stationary Gaussian stochastic process can be detected by examination of toe statistics 
of a sample from that process. It is not necessary actually to draw samples from a 
stochastic process in this investigation. A more convenient method is to use a comput- 
er to generate pseudorandom Gaussianiy distributed samples, and this was the method 
used to generate toe data presented in this report. 
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RATIOS OF MOMENTS 


A computer program that performs toe following operations is described to 
appendix A. 

1. Generation of a predetermined number of Gaussianiy distributed random 

numbers v 

/*■ 

2. Calculation of uaj directly 
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3. Calculation of directly 

2 2 

4. Evaluation of the ratio (3mg - m^J/Smg » which ideally should equal zero 

5. Repetition of steps 1 to 4 a predetermined number of times 

6. Plotting of a histogram of the values obtained from step 4 

< 

Results obtained by using this computer 
program are shown in figure 2. hi one 
case, only 50 numbers were used in the 
set of random numbers; in another case, 

200 numbers were used; and, in a third 
case, 500 numbers were used, hi each 
case, 1000 values were used in the histo- 
gram. As always in sampled statistics, 
the ideal ratio of zero was not achieved, it 
However, with sample sizes as large as 
500, the law of large numbers (ref. 3) 
applies, and the values are firmly dus- ’--if] 
tered around zero with a relatively small 
variance. These data formed the basis 
for a determination of the level of detect- 
able distortion (discussed subsequently), 
because they represent the ratios of mo- 
ments to be expected in an undistorted 
Gaussian process. 



Figure 2. - Sampled standard normal 
moments for the ratio 

(3m 2 2 - m 4 )/3m 2 2 . 
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In reference 4, information is provided from which the statistics of equation (31) 
can be calculated. 
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Thus, for very large a (n > 100), the mean of 1 - m 4 /3m 2 is zero, and its variance 
approaches 1.7/n. * 


DISTORTION 


Electrical engineers seem to be divided into two groups. One group is concerned 
with analog signals that have continuous properties, and this group of engineers speaks 
in terms of linearity, distortion, and harmonics. The other group is concerned with 
probability theory and speaks in terms of samples, distributions, and statistics. It is 
difficult to make a correlation between the two groups. Because, in this report, sam> 
pie statistics are used to Investigate linearity, a bridge between the two groups must 
be used. An attempt was made to correlate harmonic-distortion levels of sine waves 
with shifts of the statistics of Gaussian random variables. 


Distortion is defined as the percentage of content, in a signal, of the harmonics 
of a sine wave (ref. 5). If the Vj (where i * 1, 2, . . ., n) are the magnitudes of har- 
monically related sine waves in a signal, with Vj being the magnitude of the fundamen- 
tal signal, the percentage of total harmonic distortion D is 



(43) 


Distortion Is produced when a sine c 
wave is subjected to a nonlinear transfer 
and is usually an undesirable effect that - 
indicates unsatisfactory operation. In 
figure 3, the droop in the v^j/v ^ trans- 
fer curve causes the peaks of the input 

sine wave to be flattened. If the input Figure 3.- Transfer-function curve, 
were a Gaussian stochastic process of 
approximately the same size as the sine - 

wave, the input would also experience flattening of the larger values. For the ratio 
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obtained. The degree of shift would depend on the amount of distortion. Computer 
programs written to generate data for correlation of distortion levels with moment- 
ratio shifts are described in appendixes B and C. 


CORRELATION OF SIGNAL MAGNITUDES 


Processes "of approximately the same size” were referred to previously; how- 
ever, this terminology is not very precise. For the purposes of this report, the 
root mean square (rms)value of a process is used to indicate its "size. " For a sine 
wave, the rms value is 0. 707 times the peak value. For nonsinusoidal processes, the 
rms value must be calculated for each case. The rms value for a zero-mean Gaussian 
process is the square root of m g and is called the standard deviation. Obviously, the 

form factors of sine waves and Gaussian processes (i. e. , the ratios of the peak to the 
rms values) are different. However, the rms value of any given process represents the 
same amount of power as the same rms value of any other process. Furthermore, the 
rms value is linear in a linear transfer; that is, multiplication of a sine wave by a con- 
stant has the same effect on the sine-wave rms value as multiplication of a Gaussian 
process by the same constant has on the standard deviation of the Gaussian process. 


DISTORTION COMPARED WITH TRANSFER FUNCTION 


A computer program that calculates the distortion of a sine wave for a given 
transfer function is described in appendix B. Ibis program performs the following 
operations. 

1. Inputs data that specify the transfer function 

2. Fits the least-squares curve to the transfer-function data points 

■ . * 

. j . 

3. Applies a sine wave of a given rms value to the transfer function and obtains 
the output 


4. Constructs a Fourier series on the transferred signal to obtain the harmonics 

5. Calculates the distortion 

6. Repeats steps 3 to S for a different rms value, if desired 

Typical transfer functions and distortions obtained for various rms values are shown in 
figures 4, 5, and 6. 
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(c) Moment ratios. 


- r *■* .-*>*'•* . i . > 

Figure 4. - Transfer function 1 (v^ « 1. 83347 - 0. 0232 ♦ 0. 0, * to £ 7; 

v out * 7 » T in > 7), distortion for transfer function 1, and moment ratios for 
transfer function 1. 
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(a) Transfer function, -s 


(b) Distortion. 
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TRANSFERRED GAUSSIAN PROCESSES 


• A computer program that applies samples from a Gaussian process to transfer 
functions and examines the moment- ratio shifts is presented in appendix C. This com- 
puter program performs die following operations. 
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1 . 

2 . 

3 . 


Inputs the transfer-function polynomial coefficients 
Generates the random numbers 


Adjusts the sta ndar d deviation of the random numbers to some desired value 
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4. Applies the random numbers to the transfer function 

5. Calculates the first four moments of the transferred process 

2 2 

6. Evaluates the ratio (3m 2 - m^)/3m 0 

7. Repeats as many iterations as desired 

8. Plots a histogram of the results of step 6 

Histograms produced by the computer program and contrasted with the results of the 
untransferred Gaussian process are shown in figures 4(c), 5(c), and 8(c). 

DISCUSSION OF RESULTS 


A typical drooping transfer function, such as might be obtained from a "tired” 
electronics box, is shown in figure 4(a). This transfer function causes compression of 

2 2 

the input-signal peaks and should cause a positive shift in the ratio (3m 2 - m 4 )/3m 2 . 

As shown in figure 4(c), such a shift does occur. For relatively small numbers of sam- 
ples per set, the dispersion of points is so great that small amounts of distortion could 
not be detected reliably. However, for sample sets as large as 500 (fig. 4(c)), small 
amounts of distortion can be detected readily. In earlier stages of this investigation, 
a Kolmogorov-Smirnov goodness-of-fit test was incorporated into the computer pro- 
gram that is described in appendix C. However, this test was eliminated when it be- 
came evident that only rarely would a set be rejected as being from a non-Gaussian 
distribution. The moment- ratio-shift test proved to be much more sensitive, especially 
because, in the goodness-of-fit test, the sample is assumed to be from a Gaussian proc- 
ess, and the computer must have an extremely good reason before it will reject a sam- 
ple. In failure detection, the opposite assumption is more desirable, because the 
penalty for taking a good unit off-line is not high. 

The transfer function shown in figure 5(a) is the type that might be obtained in an 
amplifier with too low a power- supply voltage. The transfer function flattens very 
sharply and causes a moment- ratio shift in a Gaussian signal that has a rms value that, 
in a sine wave, would cause no distortion. This difference in distortion is caused by 
the difference between the form factors of the two signals. x The Gaussian signal has no 
rigid peak value, in theory, and always exceeds the clipping value. For a signal that 
was usually Gaussian, a standard- deviation value that would not cause excessive clip- 
ping of the signal would have to be determined, and that value would be the operating 
level. Any clipping beyond that value obviously would show up quickly. 

The transfer function shown in figure 6(a) is the type obtained when push-push 
transistors are improperly biased. The transfer function chops out the center of the 
distribution. Although not as detectable as in the other forms, small amounts of dis- 
tortion caused by this type of transfer function are still readily detectable. 


[ 


la summation, the moment- ratio- shift test appears to be an effective way of 
detecting incipient failures that are reflected in nonlinear transfer functions. The use- 
fulness of this test would be a function of the failure modes and would have to be as- 
sessed for each given situation. Nonetheless, the idea of using on-line signal statistics 
in failure detection appears to have great promise and should be pursued. 


SUGGESTIONS FOR FURTHER INVESTIGATION 


Although the results presented in the preceding sections are interesting, they 
point toward the possibility of work of a much broader and more important scope. Most 
signals used in communications are not Gaussian and are far from stationary. Voice, 
which is one of the most common types of communications signals, is neither Gaussian 
nor stationary; voice signals vary in almost every way possible. However, this does 
not mean that the use of on-line signal statistics as an aid in the detection of incipient 
failures is unfeasible. Rather, this would appear to be a situation tailormade for the 
application of learning theory, which is a discipline that has received much attention in 
recent years. 

Nonlinear operations are most likely to occur inside an electronics box. Because 
the box is probably very small (particularly in a spacecraft), both the input and output 
signals would be readily available. Furthermore, all future spacecraft will probably 
have powerful general-purpose computers on board. It seems reasonable to suppose 
that, in such a situation, the statistics of the input signal to the box could be learned 
by the computer and compared with the statistics of the output signal of the box. A 
significant shift in statistics would indicate an incipient failure. 

Several aspects of learning theory (ref. 6) will probably be required in such an 
application. The technique will most certainly require unsupervised learning (ref. 7), 
because only general characteristics will be known in advance, and the signals will not 
be stationary. Because of the nonstationary signal characteristics, a form of moving- 
window technique with optimum stopping rules (ref. 8) will probably be required. Be- 
cause any given signal has a great variety of statistics, some class of sufficient 
statistics (ref. 8) must be chosen for manipulation of each signal type. 

The most promising method of study in the use of on-line statistics to detect in- 
cipient failures appears to be the application of learning theory. Specifically, the fol- 
lowing steps should be taken. 

1. Investigation of the statistical properties of various communications signals 

2. Determination of sufficient statistics, ideally those that are generally appli- 
cable in communications signals 

3. Application of nonsupervised learning techniques, probably of a moving- 
window type, in the determination of the sufficient statistics 

4. Determination of optimum stopping rules for making M good-bajd M decisions 
about on-line equipment 



CONCLUSIONS 


The information presented in this report leads to the following conclusions. 

1. The nonlinear transfer of a Gaussian signal can be detected by using on-line 
signal statistics. 

2. The moment- ratio- shift test is an effective method for the detection of very 
small distortion levels. 

3. Further investigation is warranted, specifically on the application of learning 
theory to the problem of detecting incipient failures by using on-line signal statistics. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, July 7, 1971 
908-42-07-00-72 
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APPENDIX A 

COMPUTER PROGRAM FOR THE DETERMINATION 
OF THE RATIOS OF GAUSSIAN MOMENTS 


A computer program (Bg. A-l) that generates sets of Ggasslaaty asfarflnted 
rrv^r onAwn numbers, calculates nwxneuts of the sete, te 

moments is described in this appendix. A block diagram of the compiler program Is 
shorn in figure A-l- The operation is described as follons. 

1. The eoOtrol card {first data card) contains the following items. 

a* KNUM — ttw» muDcber of random numbers set 

^ NBAKD a w odd to initialize the random-somber generator 

c. KRDN — the number of trials to be run 

d. it — the highest order moment to be calculated 

2, The instmctica Z = EAHDOM ( KRAHD ) laiti a liae s tie 
ZOR 


3 _ The DO 6 loop determines the number of trials to be 
4. The DO 3 loop genera tes the KMPM 
«y eaaned so ta atB e calculate s 


_ 2 . rt _ 1 

«. The remainder of tte loop calculates STOR1 = (3n»2 - . 

T. TTr nnri HISI ml rpla te a and prints a histogram *be NftOH 


ntet at STOR1 


,]L 


• J \ ^ • ■ 

1 s.. 




, :.A&. 


- . ‘ 'i ‘ ’ i , - 

• - * * - >r-* T 


;Ti. 


*ynw**. 

,-W.V •••'.•fcVS •, f.pC •• ; -> y * • 

>„ ' * * "s ’ , » * - ' 





t 

c 


tax?* 

4 . 81 X 3 >* 

mu K 

ecu* ** 

Ctli* •* 

CCl^ Jl* 

can !»•* 

;ci2» jv« 

Kll- X>* 

C3xJi j** 

Kli/ X** 


V** ,aK ** B| ,na ** , »** 8 “•> 
_ w *'* •**«** *^ •» 

C v,W, 

f IS’** 

C 5 * i.'l'Ua m 
t: J 


>']"MIX< 
3 caTj-j,* 


• C^Tl'v* 

nu. Tj-nir;n> 


J= «sa;j;S5' , "Si! 

ft" - » : i; * 

& «, ■ Tr "* ! 

IZ "S 


W fc> ***** * ,c * '»”*• » c b»n«in- f •B|«V0tt|C* T*n~** , i* 






APPENDIX B 

COMPUTER PROGRAM FOR CALCULATION OF DISTORTION 


A computer program Big. B-l) that calculates fee distortion of a f^fJKwe 
applied to an input transfer function is described in this appendix. The following list 
is a descri ptio n of fee key data cards feat are used in fee operation of this computer 

program. 

1. The following control parameters are input on fee ini tial data cards. 

a. n —The number of ordered pairs (X, Y) in a set that defines fee transfer 

function 


b. NP The number of increments into which 2v is to be divided for fee 

Fourier analysis 

'.t: 

' t. NRMS — The number of nns values of the sine waves to be run 

d. bc — The order of fee polynomial to be fit to fee transfer function 

e . k — The order of toe polynomial to be fit to fee transfer function 


f. NH — The number of harmonics to be used in the Fourier series 


g. KPOL — The test to see whether data or polynomial coefficients determine 
the transfer fun c tion 

0) A sero Indicates the ordered pair (X, Y) to be read in. 

0) A one Indicates fee coefficients to be read in. 

k. CHIT — a critical value Of such esists) above which toe transfer fimction 
trices on a well-def ined value or approaches an asymptote 

L ASSXfPT — A possible asymptotic value of the transfer function 

JL Polynomial coefficient* Of each exist) are to be read in, seven to a card. 

c. 

3. Ordered pairs (g,Y) tori define a transfer function are to be read to, one 
pair per card. 

4. Therms values to be need are rend in, one per card. 

A block diagram of this computer program la shown tn figure B-l. The compwter 
p ro g ram performs the following operations. 


1. A decision is made on whether <bfa or input polynomial coefficients wUl de- 
fine toe transfer function, and fee choice to read in., - A sw ’ ti 


24 
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A block diagram of this computer program is shown in figure B-2. The computer 
program performs the following operations. 

t. A decision is made on whether data or input polynomial coefficients will de- 
fine the transfer function, and the choice is read in. 

2. The calling instructions CALL ORTHLS and CALL COEFS call library rou- 
tines that fit a K-order polynomial to the input data, if such data have been read in. 

3. The polynomial coefficients are written. 

4. The ordinates of the input data to the transfer function (if such have been read 
in) and the transferred results are written out. 

5. The DO 10 loop sets the values of the abscissas where the Fourier analysis 
will take place. 

6. The DO 20 loop sets the loop that allows iterative operation, with the rms 
value of a sine wave as the variable. 

7. The DO 30 loop (a) sets the peak value of the sine wave at 1. 414 times the 
input rms value and (b) transfers these values through the transfer function (DO 40 
loop). 


8. The calling instruction CALL DFSHIE calls a library routine that performs 
a Fourier analysis on the transferred sine wave. 


9. The remainder of the computer program calculates the sum c t the squares 
of the nonfundamental terms of the Fourier series and calculates the harmonic 
distortion. 
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Figure B-l. - Computer program lor caleaUttoo of dJetortioo. 
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APPENDIX C 

COMPUTER PROGRAM FOR CALCULATION 
OF DISTORTED GAUSSIAN MOMENTS 


A computer program (fig. C-i) that generates Gaussianly distributed pseudoran- 
dom n um bers, applies the numbers to a transfer function, and examines the ratios of 
moments of the result is described in this appendix. The following is a list of the key 
data cards used in the operation of this computer program. 

1. Initial data card 

a. NPOL — The number of polynomial coefficients to be read in 

b. NNUM — The number of pseudorandom numbers per set 

c. NRAND — An odd number that is used to initialize the random-number 
generator 

d. NRUN — The number of sets of numbers to be exercised 

e. M — The highest order moment used 

f. NRMS — The number of different standard deviations to be used 

g. INIT — A flexibility number that allows different groups of random num- 
bers to be used from a set 


h. RtCR — The increment value for a DO loop 

i. CRIT — The critical value beyond which the transfer function is well 




defined 

2. POL values — five per card , 

i. 

3. The rms values — five per card 

' H 

4. UTL1 — title for the histogram * 


$ 


A block diagram of this computer program is shown in figure C-2. The program 
operates as follows. J( 


• i 1 


t. The initial data card is read. 

- . ■ - -l 

2. The number generator is initialized by Z ■ RANDOM (NRAND). 

3. The polynomial coefficients are read. 


4. The rms values are read. 
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5. The 00 6 loop sets the outer loop for the number of number sets. 

6. The DO 3 loop generates the number sets. 

7. The DO 7 loop sets the loop to adjust the standard deviation. 

8. The DO 8 loop performs the transfer of numbers. 


9. The DO 10 loop calculates the moments and ratios of moments. 



10. The subroutine HIS1 arranges and prints a histogram of the results. 
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Figure C-l.- Computer program fofcalculatton of distorted Gaussian moments. 
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